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Abstract 

According to extensive theoretical studies of the high energy limit of QCD, inelastic interactions 
are dominated by the multi-Regge final states. The appropriate gauge-invariant objects, which 
simultaneously incorporate the transverse momentum degrees of freedom, are Reggeized gluons, 
quarks and antiquarks. In the present communication we extend part on Reggeization approach 
to Drell-Yan production of massive lepton pairs. The basic ingredient is a process of Reggeized 
quark-antiquark annihilation, QQ -> 7* — >► Z + Z - , which is described by the Reggeon-Reggeon- 
photon effective vertex E^g. We calculate transverse-momentum and invariant-mass distributions 
of Drell-Yan lepton pairs measured at the CERN SPS, FNAL Tevatron and CERN LHC in the 
different ranges of energy and rapidity. We focus on angular distributions of Drell-Yan leptons in 
different kinematical ranges. The obtained results are compared with the existing data and a good 
agreement is found. The predictions for future experiments for Drell-Yan lepton pair production 
at the CERN LHC have been made. 
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* Electronic address: nefedovma@gmail.com 
"•"Electronic address: |n. nikolaev @fz-juelich.de 
^Electronic address: saleev@samsu.ru 



2 



I. INTRODUCTION 



The Drell-Yan production of lepton pairs, the principal QCD subprocess for which is 
an annihilation of quarks and antiquarks from colliding hadrons into lepton pairs [lj], is not 
only a supplement to the standard deep inelastic (DIS) studies. Drell-Yan production stands 
out as an important source of the experimental data on the flavor content of the nucleon 
sea (see [2] and references therein). More recently there has been much discussion of the 
Drell-Yan production as a unique source of the direct information on fine spin properties 
of the proton, including transversity and a new family of transverse-momentum dependent 
(TMD) structure functions Q, 4]. 

On the discovery frontier, Drell-Yan process emerges as a background to heavy vector- 
meson production in the Standard Model (SM) and beyond, as Drell-Yan lepton pairs with 
mass above the Z— bosom mass are a major background in searches for new heavy vector 
mesons. Finally, Drell-Yan production at the LHC would probe parton distribution functions 
(PDFs) at very small x, down to x ~ 10 -6 . The experimental studies of Drell-Yan lepton pair 
production include measurements of rapidity (y), invariant mass (Q), transverse momentum 
(Qt = \qt\) of lepton pairs, and angular distributions of leptons in the virtual photon rest 
frame. 

Theoretical study of Drell-Yan lepton pair production in the perturbative Quantum Chro- 
modynamics (pQCD) framework is in a very advanced state and extends from the leading- 
order (LO) to next-to-leading-order (NLO) approximations in the strong coupling constant 
&s (5- Y] to the soft initial-state gluon resummation procedure to all orders in a s (see jsl, [9] 
and references therein). 

Ever since the pioneering works by Dokshitser et al. |_10[] and Altarelli et al. there 
has been much work of the q T — distribution of lepton pairs in the framework of standard 

in 

collinear factorization with on-mass shell partons (see [5|-l9|] and references therein). More 
recent activity on q T — distributions, which goes beyond the collinear approximation, focused 
on TMD's [121 ]. The often discussed fc^— factorization with off-shell mass partons is a part of 
such approach. At small-x of our interest the principal ingredients are off-shell properties o 
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t— channel exchanges and unpolarized unintegrated parton distribution functions (PDF) 
I4I ]. Indeed, at a deeper level of high energy pQCD, quarks, antiquarks and gluons are 
known to reggeize (l5l-ll8|. Furthermore, as has been shown by Lipatov and collaborators, 
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Reggeized gluons and quarks are the appropriate gauge-invariant degrees of freedom of high- 
energy pQCD. In the practical applications, the use of Reggeized t— channel exchanges is 
justified by the dominance of the so-called multi-Regge final states in inelastic collisions of 
high energy hadrons. 

Earlier, the quark Reggeization hypothesis has been used successfully for a description 



_Tipt 
1SK 



of different spectra of prompt photons at the Fermilab Tevatron and CERN LHC |19l-l2l|. 
electron deep inelastic scattering and prompt photon production cross sections at the DESY 
HERA Q, Q], and forward Z— boson production cross section at the CERN LHC In 
the present communication we discuss the experimental data for proton-proton and proton- 



antiproton collisions including the Fermilab Tevatron experiments 



SPS experiments 
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27|, and new data from the CERN LHC 
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[2 . 



2|j24 



25| and the CERN 



This paper is organized as follows: In Sec. II, the general formalism for description of 
Drell-Yan lepton pair production is presented. In Sec. Ill, we calculate helicity structure 
functions in the LO of Parton Reggeization Approach (PRA) and obtain master formulas for 
lepton pair spectra and angular coefficients. In Sec. IV, we present results of our calculations 
for invariant-mass and transverse-momentum spectra of Drell-Yan lepton pairs and report 
a comparison with relevant experimental data from the Fermilab Tevatron, CERN SPS and 
CERN LHC. We also compare our predictions with the data for angular coefficients and 
predict transverse-momentum dependence of the angular coefficients in the energy range of 
the CERN LHC, VS = 7 - 14 TeV. Section V contains our conclusions. 



II. DRELL-YAN PAIR PRODUCTION IN QCD AND PRA 

According to the perturbative QCD Parton Model the LO subprocess of Drell-Yan lepton 
pair production is an annihilation of quark and antiquark from the colliding hadrons into 
the virtual photon which decays onto lepton pair: 

q + q^7* + /". (1) 

In this formalism the quark and antiquark in an initial state have zero transverse momentum 
with respect to the hadron collision axis and massive lepton pair (or virtual photon) is 
produced with zero transverse momentum. The angular distribution of leptons in the virtual 
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photon rest frame should be 



dN 



1 + cos 2 9, 



(2) 



dcosO 

where 6 is the polar angle of the lepton relative to the collision axis. Experimentally, massive 
lepton pairs are produced with substantial transverse momentum q T . To describe Drell-Yan 
lepton pairs with nonzero transverse momentum, within the collinear Parton Model one 
invokes the NLO partonic subprocesses 2 — >> 2, which are of the first order in a s : 



q + q^g + ^^g + l + + l 
q + g^q + rf^q + l + + l 



(3) 
(4) 



where g is the Yang-Mills gluon. It is obvious, that in these subprocesses the lepton angular 
distribution will differ from the trivial form ([2j). If the colliding hadrons were polarized, their 
polarization transfer to the initial partons and angular distributions of the final leptons will 
change respectively. In framework of the NLO QCD and the collinear Parton Model, the 
Drell-Yan lepton pair production in collisions of polarized and unpolarized hadrons have 
been studied carefully, with exception of only the regions of small qT and Q {s-Jj. 

In the JtT— factorization approach {l3,[l^], the off-shell initial partons ((/*,#*) of nonzero 



transverse momenta are considered from the beginning. Instead of processes 
the relevant LO and NLO contributions are 



and 



g* + <f 7* / + + r 
q* + f -> g + 7* -> 9 + l+ + 1 ~ 

g* + g* -> q + 7* q + Z + + I 



(5) 
(6) 
(7) 



However, one must be careful with off-shell quarks and antiquarks as one may break the 
gauge invariance of relevant amplitudes and break the electromagnetic current conservation. 
The additional problem of the fcy— factorization calculations is a double counting, when 
subprocesses J5j), flBJ an d (ED are taken into account together. These difficulties can be 
solved in PR A, where the initial off-shell gluons and quarks are considered as Reggeons or 
Reggeized gluons and quarks, which interact with usual quarks and Yang-Mills gluons in 
a special way, via gauge invariant effective vertices which incorporate the initial and final 



state radiation effects on equal footing [T7. 
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211 ] , and inclusive prompt photon production |l9 



31 



. Our previous studies of inclusive jet 
22j have shown that in PRA it is sufficient 
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to consider only LO 2—^1 subprocess for a good quantitative description the experimental 
data. 

In PRA, the LO subprocess, which describes finite Drell-Yan lepton pairs, is an anni- 
hilation of Reggeized quark and Reggeized antiquark via virtual photon: 

Q(qi) + Q(Q2) -+ 7* ^ l + (h) + r (k 2 ). (8) 

The amplitude of the subprocess (jSJ) reads as follows 

M(QiQi = 4nae i V(x 2 P 2 )r^ Q (q 1 , q 2 )U(x 1 P 1 ) <g> Ulk^Vfa), (9) 

where is the electric charge of quark i (in units of electron charge), a is the electromagnetic 
constant, and r^g(gi,g 2 ) is the Fadin-Sherman effective vertex [171. l32(|. 



r?^ 2 ) = 7* - ^ - (io) 



x 2 S xiS 

It is easy to show that the amplitude (jSJ) is gauge invariant and ^Q^ili, 92) (?i + #2)^ = 0. 
Four-momenta of Reggeized quarks (antiquarks) have transverse components and they read 
gf = ^Pf + gf T , gf T = (0, <f iT , 0), g? = gf T = -g,| = -t { ± 0. 

The Drell-Yan pair production in the proton-proton and proton-antiproton high-energy 
collisions corresponds to the following processes 

p(Px) + p(p 2 ) / + (fcx) + r(fc 2 ) + x, (11) 
p(Px) + p(p 2 ) z+fa) + r (fea) + x, (12) 

where four-momenta of particles are shown in brackets, / = e, \i (electron or muon), q = q± + 

Q2 = ki + k 2 is the four-momentum of virtual photon, Q = yfq* and = Q 2 + g T 2 = X1X2S. 

Differential cross section for processes f lTTl ) and f [T2l ) have the standard form: 

da a 2 
d 4 qdCl = 32^ 4 SQ 4 ^ 

or 

rfa a 2 



L^W^ (13) 



dQ 2 dq^dydQ 64n 3 SQ 4LflvW>X ^ 
where y is the rapidity of virtual photon (or Z + /~ lepton pair), dfl = d(f>d cos 9 is the spatial 
angle of producing positive lepton in the rest frame of virtual photon, Pi = ^(1,0,0, 1), 
P 2 = -^(1, 0, 0, —1), is the total energy of colliding particles. Here 

IT = 2{k*k v 2 + k\k^) - Q 2 g^ (15) 



is the leptonic tensor, whereas 

W IW = J d 4 ^(P 1 P 2 |j M (^)jV(0)|P 1 P 2 ) (16) 

is the hadronic tensor. 

The convolution of hadronic and leptonic tensors reads as a sum of contributions of the 
so-called helicity structure functions Wt,l,a,aa R H| : 



7 2 

W T (1 + cos 2 9) +W L (1- cos 2 9) + 



dQ 2 dq 2 ,dydQ Mtt 3 SQ 2 

+ Wa sin 29 cos + Waa sin 2 9 cos 20 

After integration over the angles 9 and (j) in f fTTl ). we obtain 

da a 2 ( 16tt 



(17) 



dQ 2 dq 2 T dy 64n 3 SQ 2 \ 3 ' WtU ^ 
where Wtl = Wt + Wl/2. In the analysis of experimental data, the angular distribution of 
leptons is represented in terms of two sets of the angular coefficients: 

^ = (1 + cos 2 9) + A ( - - - cos 2 9 ) + A x sin 29 cos <f> + — sin 2 9 cos 2<f> (19) 

and 

dN 



dQ A + 3 
with the normalization 



One set consists of the coefficients 



(l + A cos 2 9 + n sin 29 cos <f> + ^ sin 2 9 cos 2^) , (20) 



dN\ ^ 16tt 

ST )<" = - r (21) 



* = *l, Al = ^ A2 = ^, (22) 

W TL W tl W tl 

and the other one is defined by 

2-3A) 2A 1 2A 2 , N 

A = — — , 11 — — — , y — - — . (23) 

2 + ,4 2 + A) 2 + ^o 

Helicity structure functions are obtained by the projection of hadronic tensor on the photon 
states with the different polarizations e^(g), A = ±1,0: 



W T 




(24) 


W L 


TTT fl* 1/ 


(25) 


W A 




(26) 


w AA 


= W^x^-x- 


(27) 



In the reference frame of the virtual photon its polarization 4-vector can be written in 
covariant form: 

4 1 = -L(^-^), ^ = Z", (28) 

where 4- vectors X, Y, Z satisfy following conditions: X 2 = Y 2 = Z 2 = — 1, q^X^ = q^Y 1 * = 
q^Z 1 * = 0. In the Collins-Soper frame |34j], in which we are working, these unit vectors are 
defined as: 



Z 1 * = 

x» = 



QtVs 

2Q 



q T Qr\fS 
Y" = e^T u Z a X/3 



(qP 2 )P? - (qPjP? 
(qP 2 )P? + (qPjt* 



(29) 

(30) 
(31) 



where 



rpv ( t_ pfl = _]_(pV (Q-Pj) u\ 

' Q 1 1 Vs [ 1 Q 2 q >' 



III. HELICITY STRUCTURE FUNCTIONS IN PRA 

To calculate components of the hadronic tensor W^ v and helicity structure functions 
Wt,l,a,aa we need to know squared modula of partonic amplitudes and relevant PDFs from 
colliding hadrons. In the processes of massive lepton pair production with large transverse 
momenta, where Q,qr » ^qcd — 0.1 — 0.2 GeV, one has a factorization of hard scatter- 
ing subprocesses at the scale \i ~ Qt and the perturbative QCD evolution of the parton 
distributions from some initial scale /x — 1 GeV to the hard scattering scale \i. In the 
collinear Parton Model hadronic and partonic cross sections are connected by the textbook 
factorization formula: 



da(pp — )> / + / X) = J dxi J dx2fq{xi ) ji 2 )f^{x2 ) lJ 2 )da{qq / + / ), 



(32) 



where /^(xi^, /J 2 ) is quark (antiquark) collinear PDFs. The PRA generalization of the same 
factorization formula, differential in the Reggeized partons virtualities, is 



mx u t u ii 2 )mx 2 , t 2 , ii 2 )da(QQ -> i+r). 



(33) 
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The unintegrated PDFs $^(x,t,/U 2 ) are related to their collinear counterparts f^(x,fi 2 ) by 
the normalization condition 

xf?(x,Li 2 )= I dt$ p q (x,t,Li 2 ), (34) 

which furnishes a correct transition from formulas in PRA to those in the collinear Parton 

Model. In our numerical analysis, we adopt the prescription proposed by Kimber, Martin, 

1 L 

and Ryskin (KMR) [35[ to obtain the unintegrated quark PDFs of the proton from the 
conventional integrated one. Factorization scale is chosen to be \i = ^Qt-> and varying 
£ between 1/2 and 2 serves as an estimate of the scale uncertainty. It was found, that 
magnitudes of angular observables (A , A 2) A, v) are much more stable under the scale 
variations than differential cross sections. In the kinematical regions under consideration, 
their variations under scale change are found to be below 20 %, and indicated as shaded 
bands in the figures with theoretical predictions. 

Partonic cross sections are connected with squared amplitudes of subprocesses in both 
the collinear Parton Model and PRA the standard way: 



Taking into account that 



(27r) 3 2/c 10 (27r) 3 2A: 2 o" 1 ' 



no 



d'h fd 3 h 1 f „ f 



we obtain 



Mq(mS) .> l+n - - g) |Mfa(e) g»/' + '- )|2 ^g. ( 3T) 

In the collinear Parton Model the well-known answer is 

167T 2 

3Q 4 



\M(q t q t -> = a 2 e 2 L^w™, (38) 



where 



= Xl x 2 (2P?P» + 2P?P$ - SiT) (39) 



is the quark tensor, in which it is supposed that q" = XiP[ . 
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On the other hand, in PRA we obtain for the squared amplitude of the subprocess 

16vr 2 
~3Q 



MQiQi -> /+/")| 2 = -^LoPJirw™*, (40) 



where the tensor of Reggeized quarks reads: 

<r = xix* [sr + 2(p?pz + pup?) {2xiX2S tl - h) + 

+ + tfi? ) + -(q%PZ + q^) + (41) 

4^-xi^S) 4(t 2 -ari^g) p/tp „i 

Note, the both tensors (|39l) and (KB satisfy gauge-invariance condition q^w RRA = (fw™ = 
0. It is obvious that — >• , if we put — ^ and — » x^P/ 1 . 

It is interesting to define the quark helicity structure functions %,l,a,aa respectively to 
the hadron helicity structure functions Wt,l, a, aa- Upon direct calculations we obtain 

PM _ „.,PM 



w™=Q\ w™= v£»= w™=0, (42) 



whereas 



W ™ = Q 2 + (glT±jkj! ; w m = & _ ^ t)2 (43) 

<^ = 0, W ^ = fe+M. (44) 

We can also calculate the partonic angular coefficients, which are defined like the hadronic 
angular coefficients f |22l) : 

< RA _ %r ~ Q2t) 2 wj RA _ n 2<g A _ (fir + gar) 2 .... 



P^A g 2+tl+t2 ' - ^PflA - - W PRA Q2 +tl+t2 



Upon averaging over the angle, 



n 



quarks, we obtain well-known Lam- Tung 



2, between transverse momenta of the Reggeized 



36| relation for the partonic angular coefficients 



< «o >^=< a 2 >0 . (46) 

This relation, as it will be shown below, breaks strongly for the hadronic angular coefficients 
Aq and A 2l when Q and qT become smaller at the fixed collision energy. 

The helicity structure functions W RRA at the fixed values of variables S,Q,qT,y can be 
presented via corresponding quark helicity functions w RRA : 

W™ A {S y Q,q T ,y) = ^-JdhJ dfa *i, /* 2 )$?(*2, * 2 , V 2 )wt RA - (47) 
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The experimental data for sets of angular coefficients Aq, A\, A 2 or A, /i, v are presented 
as averages over certain kinematical region of variables Q, V at the fixed value of collision 
energy \/~S. We define the relevant averages as follows 

pQmax PQTmax PUmax 

W re 9 9e {s) = J -l dQ rf / d y W ^ 9e (S,Q,q T ,y), (48) 

'J Qrnin QTmin ^ Vmin 

where J = (Q m ax ~ Qmin)(qTmax ~ qTmin)(ymax ~ Vmin) • Sometimes instead of a variable y 
one uses the Feynman variable xp, which is defined in the center of mass frame of colliding 
hadrons as follows 

XF = 7§ = % {eV ~ e ~ v) ' (49) 



and 

(ill = ; 

4g: 



dy = dXF ■ (50) 



In other cases one uses the pseudorapidity 



2 |g|-g 2 2 \l-cos0/ q 

where 

e 2r? - 



^^^Tl' I* =|^f> ^#Tf. (52) 

Our LO PRA results for the cross-sections should be corrected by the so-called K-factor, 
which includes high order (HO) QCD corrections to the LO diagrams. The main part of HO 
corrections arising from real gluon emission is already accounted in LO PRA. Still another 
part comes from the non-logarithmic loop corrections arising from gluon vertex corrections. 
Accordingly to Ref. (39)], this K-factor is written as follows 

K(QQ -> 7 *) = exp(CV^^7r 2 ), (53) 

where a particular scale choice jj? = Q^ 3 Q 2 ^ 3 for evaluation a s (jjL 2 ) in (|53l) has been advo- 



cated in [39|] . This phenomenological Anzatz can be used for large Q 2 Drell- Yan lepton pair 
production only. In case of quasi-real photon production, when q T 3> Q ~ Aqcd, we put 



K(QQ) = 1 in accordance with calculations for the real photon production 21]. The typi- 
cal numerical value of K-factor at the kinematical conditions under consideration is about 
1.3 - 1.8. 
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IV. RESULTS 



We start our comparison of theoretical predictions with the experimental data for the 
invariant-mass distributions of the Drell-Yan lepton pairs. In Fig. [U the predictions of the 
LO PRA are compared with the data from the R209 Collaboration at the two values of 
collision energy = 44 GeV and = 62 GeV. The resulting invariant-mass spectrum 
has been obtained by integration over the relevant virtual photon transverse momentum 
range and rapidity range \y\ < 4. In Fig. [2J the doubly differential cross section da/dQdy 
is presented as function of Q at the ^/~S =1.8 TeV, \y\ < 1 and < 200 GeV as it has 



been measured by CDF Collaboration 



24[]. In Fig. El the data from CMS Collaboration 



are presented as Q— spectrum normalized to the cross section in the Z— boson peak. The 
solid line represents our prediction, normalized on the value of cross section in the Z-boson 
region (<j(60 < Q < 120 GeV) = 973 pb), obtained theoretically in the Ref. {29)], see Table 
10 therein. 

As one would expect, the increase in energy improves an agreement between the theory 
and experiment. Of course, the Z— boson region, 60 < Q < 120 GeV, is beyond our 
consideration. An extension of PRA to the region of the Z— boson production demands a 
use of the unknown up today effective vertex T^g which describes the Reggeized quark- 
antiquark annihilation in Z— boson. 

To demonstrate an agreement between PRA and Drell-Yan pair spectra over the lon- 
gitudinal variables, we show in the Fig. [4] the differential cross section Q 3 da / 'dxpdQ as a 
function of the Feynman variable xp after integration over all q?. Curves 1 — 8 correspond 
to Q from 4.75 GeV till = 8.25 GeV with the step of 0.5 GeV. The data are from FNAL 
fixed target experiment |30| at the ^/~S = 38.8 GeV. 

The transverse-momentum dependence of Drell-Yan lepton pair production cross section 
is demonstrated in Figs. [5] and [6l The R209 Collaboration (26] has measured q?— spectrum 
at VS = 62 GeV, I y\ < 4, and 5 < Q < 8 GeV. We describe this data quite well, especially 
at small transverse momenta qp < 2 GeV, where the pure NLO collinear Parton Model 
calculations break down and where one usually invokes unknown (ad'hoc) nonperturbative 
intrinsic parton transverse momentum. The scale dependent uncertainties of our calculations 
are about 20 % and they become larger at q? < 2 GeV, up to 70 % at the q? ~ 0. However, 
the average value is still in agreement with experimental data. 
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The UA1 Collaboration 



271 ] data for the qr~ spectrum of Drell-Yan lepton pairs at \/~S 



630 GeV are presented as a differential invariant cross section, averaged over the virtual 
photon mass in the range of 2m M < Q < 2.5 GeV, and in the rapidity range \y\ < 1.7. 
Because the cross section grows steeply when Q —> and the lower boundary in Q is fixed 
at the minimal kinematical value Q m i n = 2m M , we need to take into account muon mass m M 
in our calculations. Formally, it boils down to an additional threshold factor (1 — Arri^/Q 2 ) 3 / 2 
in a formula for the qT~ spectrum. Taking into account the large experimental error bars, 
we can conclude that our calculations agree with the data too. 

In Fig. [?l we plot our predictions for the transverse-momentum spectra of Drell-Yan 
lepton pairs at the CERN LHC for energies y/S = 7 and y/S = 14 NeV in two regions of 
virtual photon masses: 5 < Q < 50 GeV, and 120 < Q < 200 GeV. The invariant-mass 
spectra have been obtained after integration over rapidity in the range of \y\ < 3. 

In the next section we compare our theoretical results obtained in the LO PRA with the 
experimental data for the angular coefficients in Drell-Yan pair production. We consider data 
from the NuSea Collaboration at the Tevatron Collider, which correspond to = 39 
GeV . We also make prediction for angular coefficients Aq,A 2 for CERN LHC at energies 
= 7 and VS = 14 TeV. 

NuSea Collaboration from Fermilab Tevatron recently has published data j^] for Drell- 
Yan lepton pair production in fixed-target experiment with hydrogen and deuterium targets 
and E p = 800 GeV proton beam (y/S = 39 GeV). The measurements have been done in the 
following kinematic domain: 4.5 < Q < 15 GeV, < qr < 4 GeV, < xp < 0.8. The results 
of measurements of angular distributions are presented in terms of angular coefficients A, z/, fi 
as functions of virtual photon transverse momentum. We find good agreement of our LO 
PRA calculations with data for v and A at all values of qT ) as it is shown in Figs. [8] and 
[91 Additionally, we predict /i = which is also in agreement with the data within the 
experimental error bars. 

n 

It is known since the works of Lam and Tung [37|] , that with allowance for parton sub- 
processes (J3j) and (|4|) in the NLO collinear parton model, one can obtain the relation for 
the angular coefficients, A ~ This relation is known to be valid at large g^, and it has 
been verified experimentally in the Z-boson resonance mass region in [38|. For the region 
of small dilepton masses it has not been verified yet. As it can be seen from the Figures [TU1 
and dU in our approach, this relation is approximately valid for energies below 1 TeV and 
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for the large-mass region. For large energies, and especially for the small-mass region (see 
Fig. HUJ), this relation is broken at the small values of qr- Using the formulas fj4Ul) . fj4~4"j) . (l4"7j) 
and definitions of angular coefficients (1221) . one can show that: 

52fdt& q (x u t)mx 2 ,t)x4t 
MS, Q\ y, qT = 0) = E /^ q{ximMx{C y + 2t y (54) 

q 

A 2 (S,Q 2 ,y,q T = 0) = 0, (55) 

where t — t\ = t 2j because of qyr = — fer if Or — 0. So, the value of the coefficient A at 
the qT = characterizes the smearing of unintegrated PDF in transverse momentum and it 
is increasing in the small x ~ Qt/v^ region as it shown at the Fig. [TOl 

V. CONCLUSIONS 

We reported a study of the Drell-Yan lepton pair production at LO in the Parton 
Reggeization Approach, including subprocess (|8j) with Reggeized quarks in the initial state. 
The Reggeization allows to account in a simple and compact form the initial and final state 
radiation effects with full allowance for finite transverse momenta of partons. Our theoretical 
predictions provide an adequate numerical description of a multitude of the experimental 
measurements of lepton pair distributions on the invariant mass (Q), lepton pair trans- 
verse momentum (q T ) and longitudinal scaling variable (#p) as well as lepton pair angular 
distributions at the SPS, Tevatron and LHC Colliders. This good description is achieved, 
without any ad-hoc adjustments of input parameters. By contrast, in the collinear Parton 
Model, such a degree of agreement calls for NLO and NNLO corrections and complemen- 
tary soft-gluon resummations and ad-hoc nonperturbative transverse momenta of partons. 
In conclusion, the Parton Raggeization approach has once again proven to be a powerful 
tool for the theoretical description of QCD processes induced by Reggeized gluons fusion as 
well as Reggeized quarks annihilation in the high-energy limit. 

Acknowledgements 

We are grateful to B. A. Kniehl and A. V. Shipilova for useful discussions. The work 
was supported in part by the Ministry for Science and Education of the Russian Federation 



14 



under Contract No. 14.B37.21.1182. The work of M. N. is supported also by the Grant of the 
Student's Stipend Program of the Dynasty Foundation. The work of V. S. was supported in 
part by the Russian Foundation for Basic Research under Grant 11-02-00769-a and by SFB 
Fellowship of Hamburg University (SFB-676). N. N. acknowledges a support by Institut f. 
Kernphysik, Forschungszentrum Juelich at the early stages of this study. 



[1] S.D. Drell and T.M. Yan, Phys. Rev. Lett. 25, 316 (1970). 

[2] NuSea Collaboration, L.Y. Zhu et al. Phys. Rev. Lett. 102, 182001 (2009). 

[3] J.P. Ralson and D.E. Soper, Nucl. Phys. B152 (1979) 109. 

[4] M. Anselmino, V. Barone, A. Drago, and N.N. Nikolaev, Phys. Lett. B594 (2004) 97. 

[5] W.J. Stirling, and M.R. Whalley. J. Phys. G 19, Dl (1993). 

[6] S. Gavin et al. Int. J. Mod. Phys. A 10, 2961 (1995). 

[7] E.L. Berger, J. Qiu, and X. Zhung. Phys. Rev. D 65, 034006 (2002). 

[8] J.C. Collins and D.E. Soper, Nucl. Phys. B 193 (1981) 381; J.C. Collins, D.E. Soper, and 
G. Sterman. Nucl. Phys. B 250, 199 (1985). 

[9] E.L. Berger, J-W. Qui, and R.A. Rodriguez-Pedraza. Phys. Rev. D 76, 074006 (2007). 
[10] Yu.L. Dokshitzer, D. Diakonov, and S.I. Troian. Phys. Lett. B 79 (1978) 269-272. 
[11] G. Altarelli, G. Parisi, and R. Petronzio. Phys. Lett. B 76 (1978) 351. 
[12] D. Boer, PJ. Mulders. Nucl.Phys. B569 (2000) 505. 
[13] A. Szczurek, and G. Slipek. Phys. Rev. D 78, 114007 (2008) . 
[14] A.V. Lipatov, M.A. Malyshev, and N.P Zotov. JBEP 1112, 117 (2011). 
[15] L. N. Lipatov, Nucl. Phys. B452, 369 (1995). 

[16] L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976) [Yad. Fiz. 23, 642 (1976)]; E. A. Kuraev, 
L. N. Lipatov, and V. S. Fadin, Sov. Phys. JETP 44, 443 (1976) [Zh. Eksp. Teor. Fiz. 71, 840 
(1976)]; Sov. Phys. JETP 45, 199 (1977) [Zh. Eksp. Teor. Fiz. 72, 377 (1977)]; I. I. Balitsky 
and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978) [Yad. Fiz. 28, 1597 (1978)]; Sov. Phys. 
JETP 63, 904 (1986) [Zh. Eksp. Teor. Fiz. 90, 1536 (1986)]. 

[17] V. S. Fadin and V. E. Sherman, JETP Lett. 23, 599 (1976); JETP 45, 861 (1977). 

[18] V. S. Fadin and L. N. Lipatov, Nucl. Phys. B406, 259 (1993); B477, 767 (1996). 

[19] V. A. Saleev, Phys. Rev. D 78, 034033 (2008). 



15 



[20] V. A. Saleev, Phys. Rev. D 80, 114016 (2009). 

[21] B. A. Kniehl, V. A. Saleev, A. V. Shipilova, E. V. Yatsenko, Phys. Rev. D84 , 074017 (2011). 
[22] V. A. Saleev, Phys. Rev. D 78, 114031 (2008). 

[23] F. Hautmann, M. Hentschinski, H. Jung. Preprint IFT-UAM-CSIC-12-39, LPN12-050, OUTP- 

12-08-P; e-Print: larXiv: 1205. 17591 [hep-ph]. 
[24] CDF Collaboration, A. A. Affolder et all. Phys. Rev. Lett. 84, 845 (2000). 
[25] DO Collaboration, B. Abbott et all. Phys. Rev. Lett. 82, 4769 (1999) . 
[26] R209 Collaboration, D. Antreassyan et al. Phys. Rev. Lett. 48, 302 (1982). 
[27] UA1 Collaboration, Alitti et al. Phys. Lett. B 275, 202 (1992) . 
[28] LHCb Collaboration, CERN-LHCb-CONF-2012-013. 
[29] CMS Collaboration, S. Chartchuan, JHEP 1110, 007 (2011). 
[30] P. L. Macgauhey et. al., Phys. Rev. D 50, 3038 (1994). 

[31] E. N. Antonov, L. N. Lipatov, E. A. Kuraev, and I. O. Cherednikov, Nucl. Phys. B721, 111 
(2005). 

[32] L. N. Lipatov and M. I. Vyazovsky, Nucl. Phys. B597, 399 (2001). 
[33] M. Boglione and S. Melis. Phys. Rev. D 84, 034038 (2011). 
[34] J.C. Collins, and CD. Soper. Phys. Rev. D 16, 2219 (1977). 

[35] M. A. Kimber, A. D. Martin, and M. G. Ryskin, Eur. Phys. J. C 12, 655 (2000); Phys. Rev. D 
63, 114027 (2001); G. Watt, A. D. Martin, and M. G. Ryskin, Eur. Phys. J. C 31, 73 (2003). 
[36] C.S. Lam and W.K. Tung. Phys. Rev. D 18, 2447 (1978). 
[37] C.S. Lam and W.K. Tung. Phys. Lett. B 80, 228 (1979) 
[38] CDF Collaboration, T. Aaltonen et al. Phys. Rev. Lett. 106 241801 (2011) 
[39] G. Watt, A. D. Martin and M. G. Ryskin, Phys. Rev. D 70, 014012 (2004); 70, 079902(E) 
(2004). 



16 



Q, GeV 



FIG. 1: Differential cross section of Drell-Yan lepton pair production as function of virtual photon 
mass Q. The data are from R209 Collaboration 



2J. The curve 1 - >/S = 62 GeV, the curve 2 



y/S = 44 GeV. 
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FIG. 2: Doubly differential cross section of Drell-Yan lepton pair production as function of virtual 
photon mass Q. The data are from CDF Collaboration |24j at y/S = 1.8 TeV. 
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FIG. 3: Differential cross section of Drell-Yan lepton pair production as function of virtual photon 
mass Q. The data are from CMS Collaboration [29f] at = 7 TeV. Solid line - our prediction 
normalized to the theoretical value of the cross-cection in the Z-boson region. 
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FIG. 4: Differential cross section Q 3 dcr / 'dxpdQ of Drell-Yan lepton pair production as function of 
xf integrated over all q?. Curves 1 — 8 correspond Q from 4.75 GeV till Q — 8.25 GeV with the 
step equal to 0.5 GeV. The data are from FNAL fixed target experiment 
GeV. 
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FIG. 6: Differential cross section of Drell-Yan lepton pair production as function of q?- The data 
are from UAl Collaboration Q at the \y\ < 1.7, 0.2 < Q < 2.5 GeV, y/S = 630 GeV. 



22 



10 



20 



30 



40 



50 



q T , GeV 



FIG. 7: Differential cross section of Drell-Yan lepton pair production as function of at the 
CERN LHC energy = 7 TeV, and \y\ < 3. Curve 1 - 5 < Q < 50 GeV, curve 2 - 120 < Q < 200 
GeV. 
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FIG. 8: Angular coefficient v as function of qx- The histogram corresponds to LO calculation in 



PRA with KMR 



35[ unintegrated PDFs. The data are from NuSea Collaboration 
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FIG. 9: Angular coefficient A as function of qT- The histogram corresponds to LO calculation in 



PRA with KMR 



351 ] unintegrated PDFs. The data are from NuSea Collaboration 
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FIG. 10: Angular coefficients Aq (solid curves) and A% (dashed curves) in the proton-proton 
collisions as functions of qr for the mass region 5 < Q < 50 GeV and different center of mass 
energies. Curves 1,6 - 14 TeV, 2,5 - 7 TeV, 3,4 - 2 TeV. 
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FIG. 11: Angular coefficients Aq (solid curves) and A2 (dashed curves) in the proton-proton col- 
lisions as functions of qr for the mass region 120 < Q < 200 GeV and different center of mass 
energies. Curves 1,6 - 14 TeV, 2,5 - 7 TeV, 3,4 - 2 TeV. 
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